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Abstract 


An axisymmetric, 1 10 N class, rocket configured with a free expansion between the 
rocket nozzle and a surrounding duct was tested in an altitude simulation facility. The 
propellants were gaseous hydrogen and gaseous oxygen and the hardware consisted of a 
heat sink type copper rocket firing through copper ducts of various diameters and lengths. 
A secondary flow of nitrogen was introduced at the blind end of the duct to mix with the 
primary rocket mass flow in the duct. This flow was in the range of 0 to 10% of the 
primary massflow and its effect on nozzle performance was measured. The random 
measurement errors on thrust and massflow were within +/- 1%. One dimensional 
equilibrium calculations were used to establish the possible theoretical performance of 
these rocket-in-a-duct nozzles. Although the scale of these tests was small, they simulated 
the relevant flow expansion physics at a modest experimental cost. Test results indicated 
that lower performance was obtained at higher free expansion area ratios and longer 
ducts, while, higher performance was obtained with the addition of secondary flow. There 
was a discernable peak in specific impulse efficiency at 4% secondary flow. The small 
scale of these tests resulted in low performance efficiencies, but prior numerical modeling 
of larger rocket-in-a-duct engines predicted performance that was comparable to that of 
optimized rocket nozzles. This remains to be proven in large-scale, rocket-in-a-duct tests. 


Introduction 

The rocket-based combined-cycle 
(RBCC) propulsion system is recognized 
as a promising technology for Earth-to- 
orbit vehicles 1 ' 2 . These propulsion 
systems have four modes of operation 
for the increasing speed regimes of the 
vehicle. They are: (1) rocket-ejector 
ramjet, (2) ramjet, (3) scramjet, and (4) 
rocket-only modes. Many of the 
advantages of RBCC engines result from 
certain synergistic benefits that would 
not occur if the rocket and airbreathing 
elements operated separately. One of 


these benefits is projected to occur in the 
rocket-only mode of operation at high 
altitude, when the rocket nozzle flow 
expands into the engine duct and onto 
the vehicle afterbody. This flowpath 
increases the area ratio of the rocket and, 
therefore, the specific impulse of that 
mode of operation. This rocket-in-a-duct 
operation must be considered when 
optimizing the engine flow path, because 
a well-designed ramjet or scramjet flow 
path does not necessarily result in a high 
efficiency nozzle for the rocket-only 
mode of operation. Specifically, there 
are losses in the duct due to the initial 
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shock formation as the rocket plume 
interacts with the duct and there is 
viscous dissipation in long ducts. 

Recently, a modeling effort 3 numerically 
investigated the effect of various 
parameters on rocket-in-a-duct nozzle 
performance. These parameters included: 
rocket exit area ratio (e), duct inlet area 
to rocket throat area ratio, base bleed or 
secondary flow, duct exit-to-inlet area 
ratio, and duct length-to-inlet diameter 
ratio (LTD). Results from this model 
showed that large free expansion areas at 
the rocket exit and long ducts decreased 
performance, while base bleed and larger 
duct-exit-to-inlet area ratios increased 
performance. Over the range of 
parameters investigated, the model 
predicted nozzle efficiencies from 77 to 
95% of theoretical, one-dimensional 
equilibrium, vacuum performance. 

This paper conducted an experimental 
assessment of these results using a 1 10 N 
class rocket firing into a duct. The small 
size of the test articles was expected to 
lead to increased viscous losses and 
lower efficiencies in the experimental 
data as compared to the analytical study. 
An understanding of the relevant physics 
of the nozzle flow was the goal of these 
tests. 

Apparatus and Test Procedures 

A sketch of the test configuration is 
shown in Figure 1 . It consisted of a small 
gaseous hydrogen/gaseous oxygen rocket 
of area ratio 8, located inside of a 
circular duct of diameter, D and length, 

L. The flow from the rocket nozzle 
freely expanded to the diameter of the 
duct and turned to follow the duct. A 
secondary flow of gaseous nitrogen was 


injected at the base of the rocket near the 
blind end of the duct. This secondary 
flow interacted with the free expansion 
to establish the pressure at the head end 
of the duct. This flow was predicted to 
favorably influence the nozzle 
performance in Reference 3. The entire 
assembly was mounted on a thrust stand 
and performance was measured as 
outlined below. 

Five different rocket-in-a-duct 
configurations were tested as given in 
Table 1. These consisted of three 
different free expansions at a constant 
length to inlet diameter ratio (L/D=5) 
and three different L/D ratios at a 
constant free expansion (e=8 to 23.5). 
Tests were 4 seconds long to establish 
steady state pressure at the blind end of 
the duct. No diverging duct tests were 
conducted in this series of tests. 

Secondary flow of gaseous nitrogen, 
mostly in the range of 0 to 10% of the 
core rocket flow, was introduced at the 
blind end of the duct. 

Test Articles 

The rocket injector used in this study 
was designed and fabricated by GenCorp • 
Aerojet Propulsion Division under 
contract to NASA Glenn Research 
Center as part of the Space Station 
Freedom low thrust, gaseous 
hydrogen/gaseous oxygen rocket 
technology program 4 . It was a 1 1 0 N 
thrust class platelet type injector with 
inward radial injecton of the propellants. 
Flame holding was achieved by flow 
around a bluff body attached to the spark 
plug tip. The tests were conducted at 
chamber pressures between 483 and 690 
kPa with a nominal mixture ratio of 4. 

The combustion chamber adjacent to the 
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injector face was protected with 50% 
fuel film cooling which was mixed into 
the core flow by a secondary combustion 
trip ring in a water cooled adapter. 

The rocket consisted of a heat sink type 
copper chamber, shown mounted in the 
baseline rocket test configuration 
without a duct in Figure 2. The chamber 
and throat diameters, given in Table 1, 
were 2.54 cm and 1.28 cm, respectively. 
The rocket had a straight conical nozzle 
with a 15° half angle to an area ratio e=8. 
Various copper ducts were fabricated 
from commercial tubing for these tests 
as shown in Figure 3. Their dimensions 
are given in Table 1 . They were larger in 
diameter than the rocket and had area 
ratios of 15.2, 23.5, and 33.4 based on 
rocket throat area. They were mounted 
on the large flange shown in Figure 2 for 
the rocket-in-a-duct tests. An assembly 
tool was used to ensure concentricity of 
the duct with the rocket during assembly. 
One of these assemblies is shown 
mounted on the injector and water- 
cooled adapter in Figure 4. 

Test Facility and Test Procedures 

Overview 

Testing was conducted in Research 
Combustion Laboratory, Cell 1 1 (RCL- 
1 1 ) at NASA GRC. RCL- 1 1 is a test 
•facility designed for altitude testing of 
low thrust gaseous hydrogen (GH2)/ 
gaseous oxygen (G02) rockets. The 
altitude chamber of this test facility is 
shown in Figure 5. The altitude chamber 
was 0.9 1-m in diameter with viewports 
for optical access. Mechanical access to 
the research hardware was provided by a 
roll-back capability on the tank. A two- 
stage air ejector system provided 
continuous suction to maintain a 1.4 kPa 


pressure in the tank. Rockets were 
mounted horizontally and fired into a 
water-cooled diffuser shown in Figure 4, 
during chamber roll-back. The distance 
between the diffuser and nozzle exit was 
manually adjusted by means of a bellows 
joint for different test configurations. 

The exhaust was cooled by a water spray 
and vented out through muffler stacks to 
the atmosphere. 

Data acquisition was provided by an 
autonomous, programmable, digital 
system. Readings were continuously 
scanned for out-of-tolerance conditions 
during testing. Real-time data reduction 
was performed for mass flowrates and 
performance parameters. Data storage 
was provided by floppy disks and, for 
selected parameters, by strip charts. 

More detailed information about the 
facility systems can be found in 
Reference 5. 

Thrust Stand 

The thrust stand is shown in Figure 6. 
The thrust stand was capable of making 
steady-state thrust measurements up to 
220 N thrust. The horizontally oriented 
thrust stand consisted of a mounting 
plate, a thrust plate, flexures, a 50-lbf 
measurement load cell, a 50-lbf 
calibration load cell, and a nitrogen- 
loaded ram for calibration. 

The rocket was attached to the thrust 
stand by the mounting plate. The 
mounting plate was fabricated from a 
low thermal conductivity material 
(bakelite) to minimize thermal variations 
that could be induced by the rocket. The 
mounting plate was supported by two 
lateral restraints or flexures. Three 
columns connected the mounting plate to 
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the thrust plate, which was also 
supported by two lateral restraints. 

Thrust was transmitted from the thrust 
plate to a load cell. The load cell was 
not directly connected to the thrust plate, 
but rather, had a smooth ball at its end, 
which sat against the thrust plate. This 
minimized the transmission of side loads 
to the load cells. A spring-loaded screw 
was used to maintain contact between 
the thrust plate and the flexure ball and 
to preload the measuring load cell. The 
measuring load cell was aligned axially 
with the thruster. Pressurized pistons 
applied a force against the thrust plate to 
prevent “chattering” of the ball and 
thrust plate during pressure transients 
(sea level-to-altitude, altitude-to-sea 
level) in the tank. The calibration load 
cell was attached to the ram, whose yoke 
pushed against the thrust plate during 
calibration. 

Efforts were made to avoid tares on the 
thrust stand from attachments to the 
rocket hardware. Propellant lines, water 
feed lines, and tubing for the pressure 
transducers were angled into the rocket 
hardware at 90° and contributed to 
flexure stiffness. The pressure 
transducer tubing was also looped. The 
ignition cable was contained in thin 
walled, reinforced Tygon tubing, which 
was attached to one of the columns on 
the thrust stand. 

In-situ thrust calibrations were 
performed with pressurized propellant 
lines and at altitude. Calibration loads 
were applied by the nitrogen-loaded ram, 
which pushed against the thrust plate, 
applying loads on the entire thrust stand. 
A total of 1 7 points were taken for each 
calibration, in roughly equally-spaced 


increments, going up to approximately 
35-lbf simulated thrust, then back down 
to zero. Typically three thrust 
calibrations were performed before and 
after a test series. These calibrations 
were combined to provide the thrust 
precision error for that test series. A 
straight-line fit of the last calibration 
before testing was used to calculate site 
thrust. Whenever the test tank was 
opened, a new set of thrust calibrations 
was performed. 

The nitrogen-loaded ram was also used 
to apply a simulated load on the thrust 
stand (approximately 25 lbf), while data 
was recorded by the digital data 
acquisition system, as it would during a 
real test. This was conducted before and 
after each test series to determine the 
data acquisition precision error in thrust 
(explained further in Appendix A). 

Data Reduction 

Ten basic parameters were measured in 
testing that were used in performance 
calculations: 

F s „ e \ the site force as measured by the 
thrust stand 

P r : chamber static pressure at the 
injector face 

P U mb : the ambient pressure in the 
altitude chamber 

Ph2.w- the inlet pressure to the hydrogen 

critical flow venturi 

Thj.ih- the inlet temperature to the 

hydrogen critical flow venturi 

P 02. in'- the inlet pressure to the oxygen 

critical flow venturi 

To 2 .m- the inlet temperature to the 

oxygen critical flow venturi 

P N2. in'- the inlet pressure to the nitrogen 

critical flow venturi 
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7k,., i,: the inlet temperature to the 
nitrogen critical flow venturi 
T, number ■ rocket chamber material 
temperature 

From these measured quantities and 
knowledge of the rocket-in-a-duct 
geometry, the performance parameters of 
interest can be determined: 

Vacuum Force, F nii = F siu . + P amh A t . xll 


have been tabulated for hydrogen, 
oxygen, and nitrogen in Reference 6. 

Characteristic Velocity, 

= KA.k, 

nh,„ 

Thrust Coefficient, 


P„,A 


Vacuum force was determined by free 
body diagram to be the site force plus the 
ambient pressure times the nozzle exit 
area. A eX j, is the exit area of the duct. 
Hydrogen Mass Flowrate, 


mm 




'"o' ^ 

o 

V ' JJ 


HI. in 


H2 


Oxygen Mass Flowrate, 


mm — 




f O’ ^ 

o 


P() l.m 

I If 

U 02 V 02. in 


Nitrogen Mass Flowrate, 


m v 2 = 


CA, cn <*>, 


( o A 

o 


P N 2M 


V 2 ^|f^2^l 


Total Mass Flowrate, 


m „„ = m H2+ m 02 + m n 2 


For the mass flowrate calculation, the 
discharge coefficient (Cd), the area of 
the critical flow' orifice or venturi (A ven ) 
and two real gas corrections - the ideal 
sonic flow function (0*) and the ratio of 
real to ideal sonic flow functions 
(0*/0*) were required. 0,*and 0*/0i* are 
corrections for real gas effects and are 
derived in Reference 6. 0; is strictly a 
function of the gas properties, while 
0* / 0* is also a function of inlet 
temperature and pressure. These values 


Specific Impulse, 
F. 


=■ 


m, 


Note that 



Total chamber pressure, P u „, was used in 
performance calculations. It was related 
to P ( , the static pressure measurement 
made in the injector, behind the 
combustion zone. The P, measurement 
was corrected for pressure drop across 
the combustion zone (due to momentum 
loss) and then this quantity was 
converted into stagnation or total 
pressure using the estimated velocity of 
the gases in the combustion chamber. 

The pressure drop across the combustion 
zone was measured previously in another 
copper heat sink chamber of the same 
contour as the one tested in this study. 

The rocket throat area, A h was corrected 
for thermal growth during firing, based 
on the throat temperature, T c i uimher , 
(measured by a thermocouple) and 
coefficient of thermal expansion of the 
chamber material, i.e., copper. The 
constant g c is the gravitational constant. 


NASA/TM— 1 999-209440 


5 



Uncertainty Analysis 
All measurements have an associated 
experimental error, where the error is 
defined as the difference between the 
measurement and the true value. These 
errors can be combined to provide the 
measurement uncertainty for a particular 
parameter, that is, the maximum error 
that might reasonably be expected. The 
uncertainty model for rocket testing 
recommended in Reference 7 was 
adopted for this study. The uncertainty 
analysis methodology is described in 
Appendix A. Any small bias errors were 
ignored as negligible in this analysis. 

For this study, then, no bias errors were 
assumed in the uncertainty analysis. The 
analysis is of random errors only. 

The random uncertainty of the measured 
and calculated parameters was 
determined for every test firing 
conducted in this study. There was little 
variation in the uncertainty values, from 
test to test. Typical measurement 
uncertainty values are listed in Table 2. 

The measurement uncertainties were 
generally around 1 % or less for all of the 
parameters of interest. A low 
uncertainty in vacuum thrust was 
achieved in particular, due to the fact 
that careful attention was paid to the 
thrust measurement. There were 
generally six thrust calibrations and two 
thrust data acquisition error tests 
conducted for every test series, which 
increased the sample size (degrees of 
freedom) and reduced the measurement 
uncertainty. 


Results and Discussion 

The experimental rocket performance 
parameters defined above, i.e., specific 
impulse (I sp ), thrust coefficient (Cf), and 
characteristic velocity (c*) were 
normalized to their theoretical one- 
dimensional equilibrium 8 values at the 
measured rocket chamber total pressure, 
P U)h and inlet flow composition, 
including secondary nitrogen flow. The 
one-dimensional equilibrium 8 program is 
an industry standard computer program, 
which assumes that thermodynamic 
variations occur in only one direction, 
i.e., parallel to the rocket axis. 
Equilibrium, at each axial position, is 
that composition such that the free 
energy is minimized. This composition 
results in the theoretical maximum 
performance for the given rocket 
chamber total pressure and reference 
enthalpy of the inlet flow composition. 
The normalized performance parameters 
are defined as follows: 



where the subscripts 

vac - based on vacuum thrust 

ODE - based on one-dimensional 

equilibrium 
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Values of r| lsp , r) c *. and r)cf are given in 
Table 3 for relevant numerical cases 
computed in Reference 3. The value of 
T|| Sp was given in the Reference 3. The 
value of r| c * was assumed to be 1 .0 in 
the Reference 3, but, by the definitions 
used in this paper, the secondary flow 
must be included and r|c» was reduced to 
0.926 for 8% secondary flow. This 
secondary flow had the effect of 
lowering the experimental characteristic 
velocity, c\ of the rocket chamber, since 
it did not flow through the rocket throat 
and thereby, contribute to the chamber 
pressure, P„„. The injection of secondary 
flow, therefore, was counted as a loss of 
c* in this paper. The secondary flow 
interacted with the rocket flow as it 
flowed through the duct nozzle and 
influenced in the nozzle thrust 
coefficient. By definition, thrust 
coefficient efficiency, rjcf, was 
calculated by the following equation: 



The predicted specific impulse 
efficiency, thrust coefficient efficiency, 
and characteristic velocity efficiency, 
given in Table 3, are shown in Figure 7, 
8, and 9, respectively, as a function of 
secondary flow for three different free 
expansions at L/D=5. Note the decrease 
in rjisp and r| C f as the free expansion 
percent increases and the increase in r|i sp 
and T) C f as secondary flow increases. The 
r|c* decreases with the addition of 
secondary flow as discussed above. Note 
that the increase in r|cf and the decrease 
in r|c* with the addition of secondary 
flow had canceling effects in r|i sp , but a 
slight increase was still predicted. 

The experimental data for each test run 
for this paper is given in Appexdix B. 


Data for the baseline rocket and the five 
rocket-in-a-duct configurations was 
included. The theoretical one- 

O 

dimensional equilibrium performance 
was calculated at P„„ and at the 
experimentally measured mixture ratio, 
including gaseous nitrogen diluent. 

The measured head end pressure. Pi, etl d, 
at the blind end of the duct for each test 
run was normalized to the rocket 
chamber total pressure and plotted in 
Figure 10 as a function of secondary 
flow for the five rocket-in-a-duct 
configurations tested. The curves are 
linear regression curve fits to aid in 
viewing the data. One notes that, for 
these tests, head end pressures are in the 
range of 0.2 to 1.3% of the rocket 
chamber pressure and that head end 
pressure rises with the addition of 
secondary flow. One also notes that 
larger free expansion area ratios result in 
lower head end pressures. There was 
also no discemable effect of duct length 
on head end pressure. These 
measurements indicated that head end 
pressure was controlled by the pressure 
in the free expansion zone and not by 
downstream pressure. A test at higher 
downstream pressure was conducted to 
verify that downstream pressure was not 
communicated to the pressure at the 
head end of the duct. 

The measured specific impulse 
efficiency, T|| sp , thrust coefficient 
efficiency, per, and characteristic 
velocity efficiency, r)c** given in Tables 
B2, B3, and B6, for three different free 
expansions at L/D=5 are shown in 
Figures 11, 12, and 13, respectively, as a 
function of secondary flow. The baseline 
rocket data given in Table B 1 was 
plotted for reference. Note that there was 
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undesirable scatter in r| c * shown in 
Figure 13. Ideally, this data would 
collapse to a single line as shown in 
Figure 9 for the numerical data. This 
scatter indicated that the performance of 
the injector could be sensitive to mixture 
ratio and chamber pressure, which 
influence the gaseous propellant mass- 
flows and their injection velocities. Also, 
the integrity of the seals could be 
different for each test setup, resulting in 
day-to-day variations in r|c*. This could 
explain why the baseline rocket 
performance is higher than the other data 
in Figure 13. 

Note that, as with the numerical results, 
there was a decrease in r| C f as the free 
expansion percent increases and the 
increase in rjcr as secondary flow 
increases shown in Figure 12. The T|c* 
decreased with the addition of secondary 
flow as discussed previously and had the 
undesirable scatter discussed above. 

Note, again that the increase in r| C f and 
the decrease in r|c* with the addition of 
secondary flow had canceling effects in 
T)i sp , but there was a discernable peak in 
rjisp at 4% secondary flow shown in 
Figure 1 1 . 

The measured specific impulse 
efficiency, r|i sp , thrust coefficient 
efficiency, r) C f, and characteristic 
velocity efficiency, r| c *, given in Tables 
B2, B4, and B5, for three different 
nozzle lengths are shown in Figures 14, 
15, and 16, respectively, as a function of 
secondary flow. The baseline rocket data 
given in Table B 1 was plotted for 
reference. Note again that there was 
undesirable scatter in r| c * shown in 
Figure 1 6 and that r|c* for the 3 diameter 
long duct had a low bias. Ideally, this 
data would collapse to a single line as 


shown in Figure 9 for the numerical 
data. Again, the baseline rocket 
performance is higher than the other data 
in Figure 16. 

Note that there was a decrease in r) C f as 
the length of the duct increases and the 
increase in r|cf as secondary flow 
increases as shown in Figure 15. The r| c * 
decreased with the addition of secondary 
flow as discussed previously and had the 
undesirable scatter discussed above. 

Note, again that the increase in T)ct and 
the decrease in r|c* with the addition of 
secondary flow had canceling effects in 
T|| sp , but there was a discernable peak in 
rjisp at 4% secondary flow shown in 
Figure 14. 

Finally, note that the computational data 
of Reference 3 predicted that larger 
nozzles have higher efficiencies than the 
small-scale ducts reported in this paper. 

In fact, performance of the rocket-in-a- 
duct with secondary flow was predicted 
to be comparable to that of other 
optimized rocket nozzles. Higher 
performance of larger nozzles can be 
anticipated due to the lower viscous 
losses, but this remains to be proven in a 
large-scale, rocket-in-a-duct. 

Summary 

An understanding of the relevant physics 
of rocket-in-a-duct flows was 
experimentally tested in a 1 10 N class, 
gaseous hydrogen, gaseous oxygen, 
rocket altitude test facility. Copper heat 
sink type test apparatus was employed 
and industry standard rocket test 
procedures and uncertainty analyses 
were used. These tests simulated the 
relevant flow expansion physics and 
were economical to run. Facility thrust 
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random measurement error was within 
+/- 1%. Five different rocket-in-a-duct 
configurations were tested with a 
secondary flow of gaseous nitrogen 
introduced at the blind end of the duct. 
The experimental specific impulse 
efficiency, T)i sp , thrust coefficient 
efficiency, q C n and characteristic 
velocity efficiency, r)c* for three 
different free expansions at L/D=5 and 
three different duct lengths with a free 
expansion between e=8 to 23.5 were 
given as a function of secondary flow. 
There was undesirable scatter in r|c* 
which could be injector and/or seal 
related. There was a decrease in r)cr as 
the free expansion percent increased and 


as the length of the duct increased. An 
increase in r| C f occurred as secondary 
flow increased. The T)c* decreased with 
the addition of secondary flow by 
definition. The increase in r| C f and the 
decrease in T| c * with the addition of 
secondary flow had canceling effects in 
r)| S p, but there was a discemable peak in 
r|i S p at 4% secondary flow. Prior 
numerical modeling of larger rocket-in- 
a-duct engines predicted possible nozzle 
performance comparable to that of other 
optimized rocket nozzles, but this 
remains to be proven in a large-scale, 
rocket-in-a-duct. 
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Combustion Chamber 

Throat 

Nozzle 

Duct 

length-cm 

diameter~cm 

diameter~cm 

e 

L-cm 

D-cm 

e 

%e free 

L/D 

5.08 

2.54 

1.278 

8 

24.88 

4.976 

15.2 

47.4 

5 

5.08 

2.54 

1.278 

8 

18.55 

6.185 

23.5 

66 

3 

5.08 

2.54 

1.278 

8 

30.92 

6.185 

23.5 

66 

5 

5.08 

2.54 

1.278 

8 

43.29 



66 

7 

5.08 

2.54 

1.278 

8 



33.4 

76 

5 


Table 1 . Rocket-in-a-duct configurations tested. 


Parameter 

Measurement Uncertainty (%) 

Total Chamber Pressure 

0.57 

Vacuum Force 

0.47 

Oxygen Mass Flowrate 

1.07 

Hydrogen Mass Flowrate 

1.01 

Nitrogen Mass Flowrate 

1.00 

Total Mass Flowrate 

0.44 

Mixture Ratio 

1.43 

Characteristic Velocity 

1.04 

Thrust Coefficient 

0.73 

Specific Impulse 

0.95 


Table 2. Typical Measurement Uncertainty Values 



Table 3. Computational results for rocket-in-a-duct, L/D=5, constant area. (Reference 3) 
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Appendix A: Uncertainty Analysis Methodolo2V 


Every measurement and every parameter 
calculated from measurements have an 
associated uncertainty value. Reference 
7 has a recommended methodology 
for determining the measurement 
uncertainty in rocket testing. This 
methodology was adopted for the testing 
conducted in this study and is described 
below. 

The measured parameters in this study 
for which there is an associated error are 
summarized in Table A 1 . The sources 
of error for each parameter can be 
divided into three categories: calibration 
errors, data acquisition errors, and data 
reduction errors. Each of these 
elemental sources of errors can have two 
components, a precision (random) error 
and a bias (fixed) error associated with 
it. 

The bias error is a constant, systematic 
deviation from the true value. The 
calibration of the measuring instruments 
was used to eliminate large known 
biases. Control of the measuring process 
was used to insure that there were no 
large unknown bias errors (such as 
thermal drift). Any small bias errors 
were ignored as negligible in this 
analysis. For this study, then, no bias 
errors were assumed in the uncertainty 
analysis. The analysis is of random 
errors only. 

The precision error is the variation seen 
in repeated measurements and is 
characterized by the precision error 
index, (s). The precision index for a set 
of repeated measurements is defined by: 


N f 

_ A 

z 

X> 

- X 

1=1 1 

\ 

L 


where N = number of measurements 
made, x , = the individual measurements, 

and x = the average value of the 
individual measurements. 

The sample size or the number of 
degrees of freedom (df) associated with 
the above precision index is defined as 

df = N - 1 

The calibration precision errors in the 
measuring instruments were determined 
through offline calibration (sonic 
venturis), in-situ calibration (load cell), 
or manufacturers’ specifications 
(pressure transducers, thermocouples). 
For geometric parameters (throat and 
nozzle areas), the precision error was 
estimated as the square of the specified 
diameter tolerance. For the sonic flow 
function ratios, which were estimated 
from curve fits of table data, precision 
error was determined from the Standard 


Estimate of Error (See), 



where v, = the data point, y, r = the 
equivalent calculated value from the 
curve fit, N = number of data points, and 
C = number constants used in the curve 
fit. 

Here the number of degrees of freedom 
was defined as 
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df=N-C. 

Data acquisition errors are those 
associated with the output signal from 
the measuring instrument to the 
recording device. Sources of errors here 
include the excitation voltage, signal 
conditioning, and signal gain. In this 
study, data acquisition errors were 
determined by end-to-end calibrations 
(load cell) or from the manufacturer’s 
calibration data on the data acquisition 
cards (pressure transducers and 
thermocouples). The data reduction 
errors are those due to computer 
resolution and, for this study, were 
assumed to be zero. 

The elemental sources of error for the 
measured parameters were combined 
using the root sum square. So the 
calibration precision errors (S ct ,i) and 
data acquisition precision errors (Sj,,) 
were be combined such that for a 
measuring instrument. Mi, 

S* = Jsi 7+si 

The number of degrees of freedom was 
defined as, 

M _kk±£il 

U J XU “ r* 4 

d f,a, dfju 

Performance parameters, such as 
specific impulse, thrust coefficient, and 
characteristic velocity were calculated 
from the measured parameters. 
Parameters such as vacuum force and 
mass flowrate are determined from 
measured parameters. The calculated 
parameters used in this study are 
summarized in Table A2. 


It was necessary to propagate the errors 
in the measured parameters to calculated 
parameters. The error propagation was 
approximated with Taylor’s series 
methods. Assume that the performance 
parameter, Pj, is calculated from N 
measured parameters such that, 

Pj M, M N ) 

Then the precision error for Pj was 
calculated from the precision errors of 
the measured parameters from 



The number of degrees of freedom for Pj 
can determined from 



The partial derivative, 5 Pj/dMi, is called 
the influence coefficient and is the first 
term in a Taylor’s series expansion. So, 
the partial derivative of the performance 
parameter with respect to each 
measurement in the performance 
equation must be found for error 
propagation. Propagation must be done 
with every parameter for which an error 
might exist. That is, the error of every 
calculated parameter in Table A2 must 
be related back to the error of every 
appropriate measured parameter in Table 
Al. For example, the influence of the 
precision error in the hydrogen inlet 
pressure measurement (used for 
hydrogen mass flowrate) on the specific 
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impulse error would be determined by 
the derivative of the specific impulse 
equation with respect to the hydrogen 
inlet pressure (imbedded in the total 
mass flowrate parameter). 

The value for measurement uncertainty 
for the performance parameter Pj, used 
to express a reasonable limit for error 
was defined in Reference 7 as, 

U Pj = — 1 Bpj + Spj) 

where Spj = precision error, Bpj = bias 
error, and tg$ = Students number, the 
95th percentile point for the two-tailed 
Students “t” distribution. The 
measurement uncertainty, Upj, was a 
statistical band within which the true 
value of the measured parameter was 
expected to lie, at a 95% confidence 


level. In this analysis, the bias error is 
zero. 

The Students number was used to inflate 
the value of Upj to reduce the risk of 
underestimating S/> 7 when a small sample 
was used to calculate Spj. The hs.dfpj 
value was a function of the size of the 
sample (the number of degrees of 
freedom, dfpj) used in calculating Spj. 

For small sample sizes t^dfpj will be 
large, while for larger samples tgs.dfpj 
will be smaller. The value of 
t a , = 1 .96 for an infinite number of 

V/>.oo 

degrees of freedom and f„ 5 = 2.04 at 

30 degrees of freedom, so a value of 2.0 
was assumed for tgs is used for df = 30 or 
above. 


Measured Parameters 

Chamber Pressure 

Ambient Pressure 

Site Force 

Rocket Throat Area 

Rocket Nozzle Exit Area 

Hydrogen Venturi Discharge Coefficient 

Hydrogen Venturi Sonic Flow Function Ratio 

Hydrogen Venturi Throat Area 

Hydrogen Venturi Inlet Pressure 

Hydrogen Venturi Inlet Temperature 

Oxygen Venturi Discharge Coefficient 

Oxygen Venturi Sonic Flow Function Ratio 

Oxygen Venturi Throat Area 

Oxygen Venturi Inlet Pressure 

Oxygen Venturi Inlet Temperature 

Nitrogen Venturi Discharge Coefficient 

Nitrogen Venturi Sonic Flow Function Ratio 

Nitrogen Venturi Throat Area 

Nitrogen Venturi Inlet Pressure 

Nitrogen Venturi Inlet Temperature 


Table Al: Measured Parameters 
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Calculated Parameters 

Total Chamber Pressure 

Hydrogen Mass Flowrate 

Oxygen Mass Flowrate 

Nitrogen Mass Flowrate 

Total Mass Flowrate 

Mixture Ratio 

Vacuum Force 

Characteristic Velocity 

Thrust Coefficient 

Specific Impulse 


Table A2: Calculated Parameters 
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Appendix B. Experimental Data 


Ptot(kPa) 

462.3 

462.3 

469.2 

483.0 

496.8 
517.5 

552.0 " 

558.9 
565.8 


662.4 


MR 

%N2 

Isp(sec) 

c*(m/sec) 

Cf 

Fvac(N) 

mtot(g/s) 

Isp-eff 

c*-eff 

Cf-eff 

Phead(kPa) 

Pamb(kPa) 

3.66 

N/A 

342 

2297 

1.462 

86.6 

25.8 

0.784 

0.914 

0.858 

N/A 

1.69 

3.96 

N/A 

341 

2287 

1.464 

87.4 

26.1 

0.785 

0.919 

0.855 

N/A 

1.58 

3.87 

N/A 

342 

2301 

1.456 

88.2 

26.3 

0.786 


N/A 

1.61 

3.89 

N/A 

339 

2283 

1.457 

90.8 

27.3 

0.779 



1.61 

3.82 

N/A 

346 

2335 

1.454 

93.0 

27.4 

0.795 


N/A 

1.68 

3.81 

N/A 

348 

2350 

1.452 

96.0 

28.1 

0.799 

0.939 

0.851 

N/A 

1.75 

4.01 

N/A 

339 

2295 

1.450 

102.3 

30.7 

0.781 


N/A 

1.68 

3.97 

N/A 

346 

2345 

1.449 

104.2 

30.7 

0.796 

0.942 

0.846 

N/A 

1.77 

3.94 

N/A 

345 

2331 

1.450 

105.4 

31.2 

0.793 

0.935 

0.847 

N/A 

1.84 


N/A 

342 

2332 

1.437 

121.0 

36.1 

0.786 

0.934 

0.841 

N/A 

1.80 

4.03 

N/A 

341 

2314 

1.443 

122.8 

36.7 

0.785 

0.930 

0.842 

N/A 

1.68 


Table Bl. Test data - Baseline Rocket 




669.3 

4.00 

0.0 

669.3 

3.98 

0.9 

669.3 

3.99 

2.0 

669.3 

4.07 

2.6 

669.3 



669.3 

3.98 

10.1 


676.2 4.10 0.0 



c*(m/sec) 

Cf 

Fvac(N) 

331 

2274 

1.427 

84.5 

328 

2269 

1.419 

84.3 

335 

2210 

1.488 

88.4 

335 

2207 

1.486 

88.1 


2198 

1.490 

89.6 

333 

2152 

1.517 

91.5 

329 

2284 

1.412 

86.2 

331 

2228 

1.456 

88.8 

331 

2186 

1.483 

90.4 

320 

2062 

1.520 

92.8 

329 

2262 

1.427 

101.8 

332 

2260 

1.442 

103.4 

333 

2243 

1.458 

105.0 

332 

2201 

1.480 

106.6 

337 

2224 

1.484 

107.2 

330 

2152 

1.503 

108.6 

318 

2056 

1.518 

109.8 

335 

2286 

1.438 

123.9 

335 ' 

2258 

1.455 

124.5 

334 

2222 

1.473 

126.2 

332 

2227 

1.461 

125.3 

331 

2188 

1.485 

127.5 

317 

2085 

1.490 

127.4 

330 

2292 

1.413 

122.9 


mtot(g/s) 

Isp-eff 

c*-eff 

Cf-eff 

Phead(kPa) 

Pamb(kPa) 

26.0 

0.715 

0.915 

0.781 

1.70 

1.79 

26.2 

0.708 

0.912 

0.777 

1.70 

1.71 

26.9 

0.731 

0.898 

0.816 

2.88 

1.61 

26.8 

0.732 

0.898 

0.814 

2.75 

1.66 

27.4 

0.729 

0.893 

0.816 

3.15 

1.70 

28.0 

0.734 

0.883 

0.831 

3.55 

1.68 

26.7 

0.710 



1.70 

1.79 

27.4 

0.718 

0.902 

0.796 

2.36 

1.75 

27.9 

0.722 

0.890 

0.809 

2.88 

1.75 



0.866 

0.835 

3.95 

1.57 

31.5 

0.711 

0.910 

0.781 

1.96 

1.73 



0.788 2.49 1.77 



[WltS] 


0.796 3.15 1.70 


32 


0.897 

0.808 

3.68 

1.70 

0.905 

0.812 

3.95 

1.70 




0.822 

4.34 

1.62 

35.2 


0.832 

4.87 

1.61 

37.6 

0.723 

0.918 

0.788 

3.42 

1.72 

37.8 

0.726 

0.910 

0.798 

3.82 

1.68 

38.5 

0.728 

0.900 

0.808 

4.47 

1.63 


38.5 

0.726 

39.2 

0.731 






4.60 

1.61 

5.26 

1.55 



0.875 0.820 5.80 1.62 


37.9 

0.713 

0.923 

0.773 

2.36 

1.83 


Table B2. Test data, Rocket-in-a-Duct, 6.18 cm dia., 5 dia. Long 
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Ptot(kPa) 

462.3 

462.3 

462.3 

462.3 
469.2 

469.2 
476.1 
476.1 
476.1 

476.1 
558.9 
565.8 
565.8 

~-572.7 

572.7 

~ 572.7 

579.6 

579.6 

662.4 
662.4 
662.4 
662.4 

669.3 

676.2 
683.1 
683.1 


MR 

"3.91 

4.09 

"3.91 

3.91 

"4.09 

3.90 

"4.05 

4.14 
4.20 
4.17 

4.15 
4.20 
4.20 

4.15 
4.20 
~4~16 
"709~ 

4.16 

4.09 
4.07 

4.10 
4.12 
415 
3.99 
401 
3.96 


%N2 

0.0 

o7o 

4.3 
10.0 
0.0 

2.4 
0.0 
0.4 
0.9 
1.9 

"o7 4 " 
0.0 
0.9 
0.0 
1.9 

9.8 

2.3 

4.3 
0.0 
0.4 
0.9 

1.9 

2.5 

4.6 
0.0 
9.8 


Isp(sec) 


c*(m/sec) 


328 

318 

327 

312 

"323 

"329 

325 
328" 
324 

327 

326 

322 
324 
324 

323 
"308 
"330 

322 

323 

328 
328 

327 
326 
323 

328 
311 


2269 

2190 

2177 

"2065 

2241 

2210 

2233 

2230 

"2216 

2195 

2236 

2239 

2202 

2258 
2189 
2066 
2235 
2164 
2269 

2259 
2250 
2218 
2200 
2189 
2283 
2083 


Cf 

1.419 

"7725 

1.473 

1.481 

1.413 

"77462 

"l 425 

1.442 
1.434 
1.460 
1.430 
1.411 

1.443 
1.409 
1.447 
1.464 

1.447 
1.459 
1 .395 
1.423 
1.429 
1.446 
1.452 

1.448 
1.408 
1.463 


Fvac(N) 

84.8 

84.2 

88.0 

88.4 

"78478 "" 

87.4 

87.6 
88.0 

87.7 

89.1 
102.8 

102.3 
104.5 

103.2 
106.0 

107.4 

107.2 

108.9 

119.1 

121.1 

121.9 

123.2 

124.1 

126.1 
122.8 

127.5 


mtot(g/s) 

2§3 

2§3 

27 7 

2§3 

26.8 

277) 

27.5 
27.4 

27.6 
"27.8 

32.2 

32.3 

32.8 
" 32.4 

33.4 

35.5 

33.1 

34.5 

37.6 

37.6 

37.9 
38.4 
38.8 

39.7 

38.2 
””41 .8 


Isp-eff 

”07724 

0.703 

0.736 

0.722 

0.714 

"07734 

0.718 

0.727 

0.720 

"0.730 

0.722 

0.756 

0.720 

0.716 

0.721 

0.713 

0.737 

"0726“ 

0.713 

"0726 

0.728 

0.729 

0.729 

0.728 

0.724 

0.719 


c*-eff 

0.910 
"0884 
0.890 
0.866 
0.905 
"a896 
0.900 
"0.903 
0.902 
0.896 
0.905 
0.906 
0.895 
"0.912 
0.893 
6.871 
0.910 
0.891 
0.914 
"a 9VT 
0.910 
0.901 
0.897 
0.896 
0.917 
0.872 


Cf-eff 

0.796 

"0.796 

0.827 

0.833 

0.789 

0.821 

"0.797 

"0.804 

0.799 

0.814 

0.798 

0.787 

0.804 

0.786 

0.807 

0.820 

0.809 

0.816 

0.780 

0.796 

"0.799 

0.809 

0.812 

0.813 

0.789 

0.824 


Phead(kPa) 

275 

§32 

5^93 

§36 

275 

570 

275 

§3§ 

721 

536 

447 

372 

573 

372 

§76 

751 

§36 

738 

708 

570 

§76 

725 

791 

§37 

721 

§37 


Pamb(kPa) 

i3§" 

174 

175 

T37 
T36 
749" 
749" 
T32" 

756" 
752" 
363“ 
768*" 
760 
766" 
756 
797 
745 

T3o 

784 
792 
766" 
766" 
765" 

T3o~ 

782~ 
2.27 


Table B3. Test Data, Rocket-in-a-Duct, 4.98 cm dia., 5 dia. Long 
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%N2 

Isp(sec) 

c*(m/sec) 

Cf 

Fvac(N) 

mtot(g/s) 

Isp-eff 

c*-eff 

Cf-eff 

Phead(kPa) 

Pamb(kPa) 

462.3 

4.01 

0.0 

320 

2167 

1.449 

86.4 

27.5 

0.691 

0.873 

0.793 

1.57 

1.70 

462.3 

4.06 

0.0 

317 

2166 

1.436 

85.7 

27.5 

0.685 

0.874 

0.785 

1.70 

1.73 

462.3 

4.02 

2.5 

327 

2130 

1.508 

90.1 

28.0 

0.715 

0.867 

0.825 

3.15 

1.61 

462.3 

4.10 

9.6 

311 

1962 

1.554 

92.8 

30.4 

0.704 

0.826 

0.852 

3.68 

1.70 

469.2 

3.97 

0.9 

325 

2174 

1.466 

87.7 

27.5 



0.803 

2.23 

1.64 

469.2 

4.00 

1.9 

327 

2154 

1.489 

89.0 

27.7 



0.815 

2.75 

1.60 

469.2 


■a 

323 

2085 

1.520 

91.1 

28.7 



0.832 

3.55 

1.61 

476.1 

4.03 

0.0 

323 

2207 

1.436 

87.6 

27.6 

0.698 

0.889 

0.785 

1.70 

1.75 

476.1 

4.25 

0.0 

317 

2142 

1.452 

88.3 

28.4 

0.686 

0.870 

0.789 

1.70 

1.73 

476.1 

4.09 

0.3 

322 

2172 

1.453 

88.6 

28.0 

0.697 

0.878 

0.793 

2.09 

1.67 

558.9 

4.08 

0.0 

327 

2224 

1.442 

103.6 

32.3 

0.706 

0.896 

0.788 

1.96 

2.04 

558.9 

4.09 

0.0 

325 

2198 

1.448 

103.5 

32.5 

0.702 

0.886 

0.791 

1.96 

2.02 

558.9 

4.15 

2.4 

329 

2144 

1.505 

108.1 

33.5 

0.719 

0.875 

0.822 

3.95 

1.89 

558.9 

4.12 

10.0 

311 

1996 

1.529 

109.3 

35.8 

0.705 

0.841 

0.838 

4.47 

2.03 

565.8 

4.12 

0.4 

327 

2206 

1.452 

105.4 

32.9 

0.708 

0.892 

0.793 

2.49 

2.00 

565.8 

4.23 

0.9 

325 

2180 

1.461 

106.0 

33.2 

0.706 

0.887 

0.795 

2.75 

1.95 

565.8 

4.24 

1.9 

327 

2157 

1.487 

108.0 

33.7 

0.714 

0.881 

0.810 

3.42 

1.93 

565.8 

4.36 

■a 

322 

2093 

1.511 

110.1 

34,8 

0.712 

0.868 

0.821 

4.21 

1.90 

565.8 

4.30 

KB 

325 

2095 

1.519 

110.3 

34.6 

0.718 

0.867 

0.827 

4.34 

1.90 

669.3 

4.11 

0.0 

327 

2245 

1.427 

122.6 

38.2 

0.706 

0.905 

0.780 

2.36 

2.35 

669.3 

4.05 

0.0 

327 

2237 

1.433 

122.7 

38.2 

0.706 

0.900 

0.784 

2.36 

2.33 

669.3 

4.17 

0.4 

327 

2211 

1.452 

124.8 

38.8 

0.708 

0.895 

0.792 

3.02 

2.33 

669.3 

4.11 

0.9 

331 

2224 

1.461 

125.5 

38.6 


II 1 

3.42 

2.25 

669.3 

4.10 

1.9 

332 

2193 

1.486 

127.6 

39.1 

0.723 

0.891 

0.813 

4.21 

2.24 

669.3 

4.13 

2.5 

329 

2156 

1.494 

128.4 

39.8 



4.60 

2.21 

669.3 

4.10 

4.7 

327 

2125 

1,510 

129.8 

40.4 

0.722 

0.873 

0.827 

5.13 

2.21 

669.3 

4.18 

9.8 

314 

2024 

1.521 

130.5 

42.4 

0.711 

0.853 

0.833 

5.40 

2.45 


Table B4. Test Data, Rocket-in-a-Duct, 6.18 cm dia., 3 dia. Long 
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Di 

im 

%N2 

isp(sec) 

c*(m/sec) 

Cf 

Fvac(N) 





Phead(kPa) 

Pamb(kPa) 

476.1 

3.93 

0.0 

331 

2306 

1.408 

86.5 

26.6 

0.714 

0.925 

0.772 

1.70 

1.91 

■Ban 

4.21 

0.0 

317 

2191 

1.420 

86.9 


0.686 

0.888 

0.773 

1.70 

1.75 

476.1 



323 

2252 

1.409 

86.3 

27.2 

0.698 

0.907 

0.770 

1.70 

1.86 

476.1 


0.0 

315 

2178 

1.417 

86.5 

28.0 

0.682 

0.884 

0.771 

1.70 

1.73 

476.1 


■H 

333 

2217 

1.474 

90.7 

27.7 

0.727 

0.901 

0.808 

3.28 

1.65 

476.1 

3.98 

■a 

328 

2161 

1.487 

90.9 

28.3 

0.723 

0.885 

0.816 

3.55 

1.56 

483.0 

4.09 

0.3 

325 

2234 

1.427 

88.0 

27.6 

0.703 

0.903 

0.779 

2.09 

1.70 

483.0 

4.10 

0.8 

326 

2216 

1.441 

89.4 



0.898 

0.787 

2.36 

1.75 

483.0 

4.10 

1.9 

326 

2188 

1.462 

90.4 

28.3 

0.711 

0.891 

0.798 

2.88 

1.64 


4.00 

9.8 

311 

2062 

1.480 

91.1 

29.8 

0.704 

0.866 

0.813 

3.82 

1.58 

565.8 

4.24 

0.0 

322 

' 2247 

1.405 



0.697 

0.911 

0.765 

2.09 

1.83 

572.7 



326 

' 2240 

1.427 

104.5 

32.7 

0.707 

0.909 

0.777 

2.49 

1.70 

572.7 

4.25 

0.9 

324 

2221 

1.432 

105.5 

33.2 

0.704 

0.904 

0.779 

3.02 

1.66 

572.7 

4.34 

1.9 

322 

2180 

1.450 

106.9 

33.8 

0.703 

0.894 

0.788 

3.55 

1.61 

579.6 

4.17 

0.0 

324 

2266 

1.401 

103.6 

32.6 

0.700 

0.916 

0.764 

2.09 

1.75 

579.6 

4.06 

2.5 

329 

2221 

1.454 

108.0 

33.4 



0.796 

3.95 

1.50 

579.6 

4.29 

mm 

321 

2139 

1.472 

108.9 



0.885 

0.802 

4.34 

1.58 

579.6 

4.31 

EH 

305 

2032 

1.471 

108.9 

36.4 

0.691 



4.60 

1.66 

669.3 

3.98 

0.0 

325 

2283 

1.396 


37.5 

0.701 

0.916 

0.766 

2.49 

1.54 

669.3 

4.11 

0.8 

323 

2231 

1.420 


38.6 

0.700 

0.902 

0.776 

3.42 

1.52 

669.3 



324 

2201 

1.443 

124.3 

39.1 

0.706 

0.895 

0.789 

4.21 

1.55 

676.2 

4.05 

0.0 

324 

2272 

1.397 

121.6 

38.3 

0.699 

0.914 

0.765 

2.49 

1.45 



0.4 

324 

2260 

1 .404 

121.4 

38.2 

0.701 

0.911 

0.768 

3.02 

1.50 

676.2 

4.14 

0.9 

321 

2224 

1.417 

122.4 

38.8 

0.696 

0.901 

0.774 

3.42 

1.48 

676.2 

4.05 

2.6 

328 

2224 

1.447 

124.9 

38.8 



0.793 

4.60 

1.45 

676.2 

4.12 

■HI 

321 

2167 

1.454 

126.5 

40.1 



0.796 

5.00 

1.59 

676.2 

4.12 

9.8 

307 

2062 

1.461 

126.7 

42.0 

0.695 

0.868 

0.801 

5.53 

1.69 


Table B5. Test data, Rocket-in-a-duct, 6.18 cm dia., 7 dia. long 
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Ptot(kPa) 

476.1 

476.1 
483.0 
483.0 " 
483.0 
483.0 

483.0 
489^9 
489.9 
489.9 
489.9 
572.7 
579.6 
579~6 
579.6 
579.6 
586.5 
586.5 
586.5 
669.3 

676.2 
676^2 
676.2 
6834 

683.1 
683.1 
683.1 


MR 

4.16 

4735 " 

3.88 

4.22 

4.09 

4.04 

4.04 

4.04 

4.06 

4.08 

TicT 

4.09 
4.26 

4.13 

4.14 
4.20 
4.32 
4.32 
4.31 
4.04 

4.08 
4.11 

4.09 
4.02 

4.11 

4.12 

4.15 


% N 2 

0.0 " 

0.0 

O.O" 

0.0 

0.4 

0.9 

1.9 

2.3 " 

4.2 
9.5 
14.0 
0.0 
0.0 
0 . 4 _ 

“ 09 " 

T 9 " 

2.3 " 

4.3 

9.4 
0 . 9 ~ 
0.0 
19 " 

2.5 
0.0 
0 A~ 
4 . 6 "" 

TeT 


Isp(sec) 


c*(m/sec) 


319 

" 3 T 5 

"330 

*328 

330 

332 

327 

323 

3 U 

304 

327 

Tm 

329" 

329 " 

328 " 

324 ” 

323 

7 TT 

337 

327 

329 “ 

329 " 

332 

337 

" 325 " 

313 " 


2177 

2186 

2281 

2176 

2227 

2234 

2227 

2169 

2114 

2024 

1941 

2265 
2212 
2246 
2237 
2190 
2164 " 
2122 
2028 
2259 
2264 
2230 
2200 
2296 

2266 
2155 
2048 


Cf 

1.436 

7.415 
1.417 
1.425 
7444 
1448 
1.465 
1 .478 
1500 
1.522 
1.534 
1.422 
1.412 

1.437 
1.444 

1.468 

1.469 
1.491 
T 516 
1.439 

77T7 

1 . 447 " 
1.467 

1.416 
1 .433 
1.477 
1.497 


Fvac(N) 

87.1 

86.3 

88.1 
88.1 

89.0 

89.6 

_____ 

92.3 

94.2 

95.8 

96.7 

105.0 

105.2 
1 06/4 

107.2 
108.9 
1 10.6 

112.3 

114.4 

123.8 

123.5 

125.3 

1277 “ 

123.7 

125.5 

129.3 
131.2 


mtot(g/s) 

279 

279 

272 

28 A 

27 6 

277 
279 
2^8 
297 

371 

3 Z 5 
326 
337 
329 
33 2 
337 
377 
3^4 

372 
387 
387 
387 
397 
387 
387 
406 
427 


Isp-eff 

0.679 

0 . 672 " 

0.702 

0.673 

0.699 

0 . 705 " 

0.712 

0.703 

0 . 70 T 

0.699 

0 . 693 " 

0.698 

0.677 

0.701 

0.703 

0.704 

0.698 

0 . 702 " 

0.700 

0.707 

0 . 695 " 

0.706 

0.708 

0 . 706 " 

0.705 

0.706 

0.697 


c*-eff 

oTTf 

0.891 

0.914" 

0.882 

0.901 

0.904 

0.902 

0.883 

0.868 

0.851 

0.834 

0.913 

0.897 

0.909 

0.907 

0.894 

0.888 

0.878 

0.858 

0.912 

0.912 

0 . 906 " 

0.896 

0.922 

O . 91T 

0.886 

0.862 


Cf-eff 

0.770 

0.755 

0.767 

0.763 

0.776 

0.780 

0.789 

0.796 

0 . 809 " 

0.822 

0.830 

0 . 765 " 

0.756 

0.772 

0.776 

0.788 

0.786 

0.799 

0.814 

0.776 

0.763 

0.779 

0.791 

0.764 

0.771 

0.796 

0.808 


Phead(kPa) 

117 

130 

130 

TT 7 

144 

157 

183 

209 

223 

262 

275 

244 

144 

Tto 

196 

223 

249 

T02 

328 

223 

170 

2~62 

202 

170 

209 

042 

202 


Pamb(kPa) 

164 

279 

286 

164 
28 ) 
175 
170 
161 
155 
154 
160 
184 
175 
179 
175 
170 
161 
161 
166 
172 
174 
162 
163 
181 
181 

165 
170 


Table B6. Test data, Rocket-in-a-duct, 7.38 cm dia., 5 dia. long 
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Figure 1 . Sketch of test apparatus 


( ^Rocket^ ) 


Duct Mounting Flange 





Figure 2. Heat sink type copper rocket in baseline rocket test configuration 
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Figure 3. Various copper ducts fabricated for the test series from commercial pipe 



Figure 4. Rocket-in-a-duct test configuration 
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Figure 5. Research Combustion Laboratory Cell-1 1 



Figure 6. Horizontally-oriented thrust stand 
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Secondary Flow ~ % 


Figure 7. Predicted specific impulse efficiency of a rocket-in-a-duct as a 
function of secondary flow and free expansion area ratio for L/D=5, constant 
area duct. (CFD data from Reference 3.) 
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Secondary Flow ~ % 


Figure 8. Predicted thrust coefficient efficiency of a rocket-in-a-duct nozzle 
as a function of secondary flow and free expansion area ratio for L/D=5, 
constant area duct. (CFD data from Reference 3.) 
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Figure 9. Characteristic velocity efficiency of a rocket-in-a-duct as a 
function of secondary flow and free expansion area ratio for L/D=5, constant 
area duct. (CFD data from Reference 3.) 
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Figure 1 0. Ratio of duct head end pressure to rocket chamber total pressure. 
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Figure 1 1 . Specific impulse efficiency of rocket-in-a-duct for L/D=5 
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0.88 



Figure 12. Thrust coefficient efficiency of rocket-in-a-duct for L/D=5. 
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Secondar 


Figure 13. Characteristic velocity efficiency o 
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rocket-in-a-duct for L/D=5. 






Figure 14. Specific impulse efficiency of rocket-in-a-duct with free 
expansion £=8-23.5 as a function of percent secondary flow and length of 
duct 
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Secondary Flow ~ % 


Figure 15. Thrust coefficient efficiency of rocket-in-a-duct with free 
expansion £=8-23.5 as a function of percent secondary flow and length of 
duct 
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Figure 1 6. Characteristic velocity efficiency of rocket-in-a-duct with free 
expansion £=8-23.5 as a function of percent secondary flow and length of 
duct 
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